Novel features of J/ ^ dissociation in matter 
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We make a detailed study of the effect that the recently predicted modification of the in-medium masses 
of charmed mesons would have on J/\& dissociation on pion and p- meson comovers in relativistic heavy ion 
collisions. We find a substantial dependence of the J/\& absorption rates on the density and temperature of the 
nuclear matter. This suggests that a quantitative analysis of J/\t dissociation in nucleus nucleus collisions should 
include the effects of the modification of meson masses in dense matter. 



The modification of hadronic properties in a 
nuclear medium may be related to the partial 
restoration of chiral symmetry Jl]] - an idea 
which is currently receiving considerable atten- 
tion. Some experimental evidence for such effects 
has been discovered recently - e.g., see Ref. || 
for a review. 

There is also a great deal of interest in pos- 
sible signals of Quark-Gluon Plasma (QGP) for- 
mation (or precursors to its formation) and J/ty 
suppression is a promising candidate which has 
recently shown an anomalous result |jpj]. On 
the other hand, there may be other mechanisms 
which produce an increase in J/vE' absorption in 
a hot, dense medium. We are particularly inter- 
ested in the rather exciting suggestion, based on 
the quark-meson coupling (QMC) model ||, that 
the charmed mesons, D, D, D* and D* , should 
suffer substantial changes in their properties in a 
nuclear medium This might be expected to 
have a considerable impact on charm production 
in heavy ion collisions. 

In Ref. § it was found, for example, that at 
a density pB=3po(po=0.15 fm -3 ) the D-meson 
would feel an attractive scalar potential of about 
120 MeV and an attractive vector potential of 
about 250 MeV. These potentials are comparable 
to those felt by a if - -meson plpl, while the total 
potential felt by the D is much larger than that 
for the vector mesons, p, lu and (f> Within 
QMC it is expected that the mass of the J/'i' 



should only be changed by a tiny amount in nu- 
clear matter @||. A similar result has also been 
obtained using QCD sum rules Jl(J. 

In the light of these results, it seems that the 
charmed mesons, together with the K~ , are prob- 
ably the best candidates to provide us with infor- 
mation on the partial restoration of chiral symme- 
try. Both open charm production and the disso- 
ciation of charmonia in matter may therefore be 
used as new ways of detecting the modification of 
particle properties in a nuclear medium. 

The suppression of J /iff production observed in 
relativistic heavy ion collisions, from p+A up to 
central S+U collisions, has been well understood 
in terms of charmonium absorption in the nuclear 
medium. However, recent data from Pb+Pb col- 
lisions show a considerably stronger J/ VP suppres- 
sion ||. In an attempt to explain this "anoma- 
lous" suppression of J/^S production, many au- 
thors have studied one of two possible mecha- 
nisms, namely hadronic processes pd]Jl2] , |h| . |l5| 
and QGP formation (l6| (see Ref. P| for a re- 
view) . 

In the hadronic dissociation scenario |ll[] it is 
well known that the J/'I' interacts with pions 
and p-mesons in matter, forming charmed mesons 
through the reactions, n+ J/^—>D* +D,D* + D 
and p+J/'ty—fD + D. The absorption of J/* 
mesons on pions and p-mesons has been found 
to be important (see Refs. [p| p^Jl7[ | and refer- 
ences therein) in general and absolutely neces- 
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sary in order to fit the data on J production. 
Furthermore, the absorption on comovers should 
certainly play a more important role in S + U 
and Pb + Pb experiments, where hot, high den- 
sity mesonic matter is expected to be achieved. 

J /if! dissociation on comovers, combined with 
the absorption on nucleons, is the main mecha- 
nism proposed as an alternative to that of Mat- 
sui and Satz |ll| - namely the dissociation of the 
J /if! in a QGP. Note that both the hadronic and 
QGP scenarios predict J/vP suppression but no 
mechanism has yet been found to separate them 
experimentally. 
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Figure 1. ir+J/if! (a) and p+J/if! (b) dissociation 
cross sections as functions of the invariant colli- 
sion energy, s 1 / 2 . Results are shown for vacuum 
(the dotted line), po (the dashed line) and 3po 
(the solid line). 



Within the hadronic scenario the crucial point 
is the required dissociation strength. In par- 



ticular, one needs a total cross section for the 
n, p+J/if! interaction of around 1.5-^3 mb in or- 
der to explain the data in heavy ion simula- 
tions . Recent calculations Jl5| of the reactions 
n+J/Hf^D + D*,D + D* and p+J/f>^D + D, 
based on D exchange, indicate a much lower cross 
section than this. 

The main uncertainty in the discussion of the 
J /if! dissociation on a meson gas is given by the 
estimates of the 7r, p+J/if! cross section |L^|. The 
predictions available for the n+J/if! cross sec- 
tion are given in Refs. Following the 
meson exchange model of Rcf. [^5|, we show by 
the dotted line in the Fig. [l] the ir+J/if! (a) 
and p+J/if! (b) dissociation cross sections cal- 
culated in free space. Moreover, the upper axis 
of Fig. 0a) indicates the 7r-meson kinetic energy, 
Tjr, given in the J /if! rest frame, which indicates 
that the pions should be sufficiently hot to be 
above the DD* production threshold and to dis- 
sociate the J/^-meson. By taking a thermal pion 
gas with average T„-~I50 MeV, one might con- 
clude that independent of the ir+J/if! dissocia- 
tion model used |p^ , |l5| , the rate of this process 
is small. However, this situation changes when 
the in-medium potentials of the charmed mesons 
are taken into account, because they lower the 
iv+J/if!—>D+D* threshold. The upper axis of 
Fig. |l|b) shows the p-meson kinetic energy T p 
in the J/ if! rest frame. As was discussed in 
Refs. Ugljl, the J/* dissociation might pro- 
ceed on thermal p-mesons, because of the low 
p+J/^^D+D* reaction threshold. 

As far as the meson properties in free space are 
concerned, the Bethe-Salpeter (BS) and Dyson- 
Schwinger (DS) approaches have been widely 
used p8[ . The application of BS approach to the 
description of heavy-light quark systems allows 
one to describe the D and B meson properties in 
free space quite well The DS approach at 

finite baryon density was used |20j for the calcu- 
lation of the in-medium properties of p, u) and <f> 
mesons. The modification of the p and to meson 
masses resulting from the DS equation is close to 
the calculations with the quark-meson coupling 
(QMC) model |,||[ . while the <^-meson mass re- 
duction from Ref. 20| is larger than the QMC 
result. 
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Here, we use the quark-meson coupling 
model ||, which has been successfully applied 
not only to the problems of nuclear binding and 
charge densities but also to the study of 
meson properties in a nuclear medium 
A detailed description of the Lagrangian den- 
sity and the mean-field equations are given in 
Refs. §,@J|J|l|. The Dirac equations_ for the 
quarks and antiquarks in the D and D meson 
bags (q stands for the light quarks hereafter), ne- 
glecting the Coulomb force, are given by 0: 



tpu(x) 



= 0, 



H-d x - (m, - V q ) T 7° (Vj - \v^j 



= 0, 



[ij ■ d x - m c ] ip c (x) (or ip s (x)) = 0. 



(1) 



(2) 
(3) 



The mean-field potentials for a bag in symmetric 
nuclear matter are defined by V q =g q a, V q = g^cu 
and V q =g q p b, with g q , g^ and g q p the correspond- 
ing quark and meson- field coupling constants. 

The normalized, static solution for the ground 
state quarks or antiquarks in the meson bags may 
be written as 0: 

i) f {x)=N f e- i *f t ' R H f {x), (j = D,D), (4) 

where f—u, u, d, d, c, c refers to quark flavors, 
and Nf and ipf(x) are the normalization factor 
and corresponding spin and spatial part of the 
wave function. The bag radius in medium, R*, 
which generally depends on the hadron species to 
which the quarks and antiquarks belong, is de- 
termined through the stability condition for the 
mass of the meson against the variation of the bag 
radius . The eigenenergies e/ in Eq. (Q) in 

units of are given by 



n* q ±R*[v« + -v« 



= 0* ± R* | I 

— 6^ — Qr.t 



2 p 



(5) 

(6) 
(7) 



where 0*=J x 2 +(R*m*) 2 , with 



^</-fl , |c^ 



and i} c —^Jx 2 +(R*m c ) 2 . The bag eigenfrequen- 

cies, x q and determined by the usual, lin- 

ear boundary condition . 

The D and D meson masses in symmetric nu- 
clear matter are given by: 

= -^—f^ (8) 



<9m 



D 



0. 



(9) 



In Eq. (g), the Z£> parametrize the sum of the 
center-of-mass and gluon fluctuation effects, and 
are assumed to be independent of density plj . 
The parameters are determined in free space to 
reproduce their physical masses. 

In this study we chose the values 
m q =m u —md=5 MeV and to c =1300 MeV for the 
current quark masses, and i?jv=0.8 fm for the bag 
radius of the nucleon in free space. Other input 
parameters and some of the quantities calculated 
are given in Refs. We stress that while the 

model has a number of parameters, only three of 
them, g q , g^ and g q , are adjusted to fit nuclear 
data - namely the saturation energy and den- 
sity of symmetric nuclear matter and the bulk 
symmetry energy. Exactly the same coupling 
constants, g%, and g q , are used for the light 
quarks in the mesons as in the nucleon. Through 
Eqs. (]l]) - (0) we self-consistently calculate ef- 
fective masses, m* D , and mean field potentials, 
V§ui pi m symmetric nuclear matter. The scalar 
and vector potentials felt by the D and D mesons 
are given by H : 



U D 



U s =i 
T(V q 



>D 
1 

" 2 



m D , 



(10) 
(11) 



where, V q — 1.4 VJ, which is assumed to be the 
same as that for the K + and K~ mesons HQ. 
The isovector meson mean field potential, V q , is 
zero in symmetric nuclear matter. 

Finally, the D, D* and p- meson potentials used 
in further calculations are shown in Fig. || as a 
function of the baryon density, in units of po=0.15 
fm -3 . Note that these potentials enter not only in 
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Figure 2. The scalar (a) and vector (b) potentials 
for the D, D* and p mesons, calculated for nu- 
clear matter as functions of the baryon density, in 
units of the saturation density of nuclear matter, 
/?o=0.15 fm -3 . Scalar potentials for D and D* 
are indistinguishable. 
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the final state phase space (which becomes larger 
since the scalar masses are reduced in matter), 
but also in the reaction amplitude and the initial 
p-meson mass. Furthermore, as observed earlier, 
the properties of the J/ty meson are not signifi- 
cantly altered in medium within QMC. 

Note that the total D~-meson potential is re- 
pulsive, while the D + potential is attractive, 
which is analogous to the case for the K + and 
K~ mesons, respectively HQ]. The threshold re- 
duction is quite large when the nuclear density 
becomes large for the D + D~ pairs. Note that a 
similar situation holds for the K + and K~ pro- 
duction and, indeed, enhanced i4T~-meson pro- 
duction in heavy ion collisions, associated with 
the reduction of the production threshold, has 
been partially confirmed experimentally P,|22| . 



We first discuss the thermally averaged cross 
sections, (av), for n+J/ty and p+J/^ dissocia- 
tion in Figs. H and [|, when the free masses are 
used for the charmed mesons and the in-medium 
potentials are set to zero. They are shown by the 
dotted lines. These results are needed for com- 
parison with the calculations including the po- 
tentials. 
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Figure 3. Thermally averaged tt+ J/\& absorption 
cross section as a function of the pion gas temper- 
ature, T, and chemical potential, /i, at pj/^=0 
(a) and as a function of the J/ ^-momentum at 
T=180 MeV (b). The results are shown using the 
same notation as in Fig. 1 . The shadowed area in- 
dicates the temperatures expected to be achieved 
in heavy ion collisions. 



Since the pions are almost in thermal equilib- 
rium, their energy spectrum is given by a Bosc 
distribution with temperature, T, and chemi- 
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cal potential, p, where we have used the value, 
p=126 MeV @. The thermally averaged tt+ J/* 
cross section can be obtained by averaging over 
the 7r-spectrum at fixed J/^-momentum. The 
dotted line in Fig. ||a) shows (av) as a function 
of the pion gas temperature, T, which was cal- 
culated with zero J/vE' momentum relative to the 
pion gas. 

The shadowed area in Fig. |]a) indicates the 
temperature range corresponding to the pion den- 
sities 0.2 — 0.8 fm -3 , which are expected to be 
achieved in the heavy ion collisions presently un- 
der consideration. In vacuum the n+J/ty dissoci- 
ation cross section is less than about 0.3 mb. The 
thermally averaged absorption cross section for 
temperature, T=180 MeV, is shown in Fig. ||b) 
(the dotted line) as a function of the J/\P momen- 
tum. The thermally averaged cross section which 
we find, (av), would be very difficult to detect 
with the present experimental capabilities. 

A similar situation holds for the p+J/^ disso- 
ciation, as illustrated by the dotted line in Fig. |[ 
Indeed, the J /\& absorption on comovers seems to 
be negligible jl5) in comparison with that needed 
to explain the J/\I / suppression, provided that 
we use the vacuum properties of the charmed 
mesons. 

On the other hand, this situation changes dra- 
matically when we consider the effect of the vec- 
tor and scalar potentials felt by the charmed 
D, D* and p mesons as calculated by Eqs. (10) 
and (10). The cross sections calculated for 



7r, p+ J/\E' collisions with the in-medium poten- 
tials are shown in Fig. [I], for densities, po (the 
dashed line) and 3po (the solid line). The dot- 
ted line in Fig. [I] indicates the free space cross 
sections. 

Clearly the J/\E' absorption cross sections are 
substantially enhanced for both the tt+J/^> and 
p+J/if> reactions, not only because of the down- 
ward shift of the reaction threshold, but also be- 
cause of the in-medium effect on the reaction am- 
plitude. Moreover, now the J/ty absorption on 
comovers becomes density dependent - a crucial 
finding given the situation in actual heavy ion col- 
lisions. These effects have never been considered 
before. The comover absorption cross section is 
calculated as a function of baryon density for the 
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Figure 4. Thermally averaged p+J/ty absorption 
cross section as a function of the p-meson gas tem- 
perature T with pj/*=0 (a) and as a function of 
the J/f-momentum at T=180 MeV (b). Nota- 
tions are similar to Fig. [| 



first time. 

The thermally averaged, in-medium w+J/ty 
and p+J/^ absorption cross sections, (av), are 
shown by the dashed and solid lines in Figs. || 
and |], respectively. We find that (av) depends 
very strongly on the nuclear density. Even for 
Pj/*=0, with a pion gas temperature of 120 MeV, 
which is close to the saturation pion density, the 
thermally averaged J /if? absorption cross section 
on the pion, at ps=3po, is about a factor of 7 
larger than that at ps=0 (i.e., with no effect of 
the in- medium modification - see Fig. [|a)). 

The thermally averaged p+J/^f? dissociation 
cross section at p_B=3po becomes larger than 
1 mb. Thus, the J/^f! absorption on p-mesons 
should be appreciable, even though the p-meson 
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density is expected to be small in heavy ion col- 
lisions. We note that dynamical calculations Q 
suggest that the p-density should be around half 
of the pion density in Pb+Pb collisions. 

In order to compare our results with the 
NA38/NA50 data |@ on J/* suppression in 
Pb+Pb collisions, we have adopted the heavy ion 
model proposed in Ref. jl2| with the Et model 
from Ref. |D| . We introduce the absorption cross 
section on comovers as function of the density of 
comovers, while the nuclear absorption cross sec- 
tion is taken as 4.5 mb |13| ], Our calculations 
are shown in Fig. |j| The dashed line in Fig. || 
shows the calculations with the phenomenological 
constant cross section for J absorption on co- 
movers (av)~l mb and is identical to the results 
given in Ref. |13|. The solid line in Fig. || shows 
the calculations with the density dependent cross 
section (av) for J/*? absorption on comovers cal- 
culated in this work. Both curves clearly repro- 
duce the data || quite well, including most recent 
results from NA50 on the ratio of J/<3/ over Drcll- 
Yan cross sections, as a function of the transverse 
energy up to Et=W0 GeV. It is important to 
note that if one neglected the in-medium modi- 
fication of the J/*f? absorption cross section the 
large cross section (eru)~l mb cannot be justified 
by microscopic theoretical calculations and thus 
the NA50 data cannot be described. 

Furthermore, we notice that our calculations 
with in-medium modified absorption provide a 
significant improvement in the understanding of 
the data || compared to the models quoted by 
NA50 [e3. The basic difference between our re- 
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24 1 is that in pre- 
12] the cross sec- 



suits and those quoted by NA50 
vious heavy ion calculations [^|4 
tion for J absorption on comovers was taken as 
a free parameter to be adjusted to the data pj-|2~| 
and was never motivated theoretically. 

To summarize, we have studied J/vP dissocia- 
tion in a gas of 7r and p mesons, taking into ac- 
count for the first time the density dependence of 
the scalar and vector potentials which the mesons 
feel in nuclear matter. We found a substantial 
density dependence of the J/ 1 ^ absorption rate 
as a result of the changes in the properties of 
the charmed mesons in-medium. This aspect has 
never been considered before when analyzing J/\& 
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Figure 5. The ratio of the J/^> over Drell-Yan 
cross sections from Pb+Pb collisions as function 
of the transverse energy Et- Data are from 
Ref. The solid line shows our calcula- 

tions with the density dependent cross section for 
J/vE' absorption on comovers. The dashed line 
indicates the calculations with phenomenological 
cross section (av)c^l mb fL3[ . For both calcu- 
lations the nuclear absorption cross section was 
taken as 4.5 mb. 



production in heavy ion collisions. 

Moreover, when we introduce the density de- 
pendent cross section on comovers into a heavy 
ion calculation and compare our results with most 
recent NA50 data |||4| , our calculations indicate 
very good agreement with the NA50 data up to 
transverse energy of order 100 GeV. 
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